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Abstract 

Recently, polarons in the Peyrard-Bishop-Holstein model under DC elec- 
tric fields were established to perform Bloch oscillations, provided the charge- 
lattice coupling is not large. In this work, we study this model when the 
charge is subjected to an applied field with both DC and AC components. 
Similarly to what happens in the rigid lattice, we find that the carrier under- 
goes a directed motion or coherent oscillations when the AC field is resonant 
or detuned with respect to the Bloch frequency, respectively. The electric 
density current and its Fourier spectrum are also studied to reveal the fre- 
quencies involved in the polaron dynamics. 
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1. Introduction 

In perfect crystals whose lattices generate a periodic potential, Bloch the- 
orem predicts uniformly states extended over the whole system |l| . However, 
if a uniform electrical field F is applied, in addition to this periodic poten- 
tial, all states become spatially localized due to Bragg reflections 0, 0, S]- 
In such a case, neglecting scattering effects, electron keep oscillating within 
a finite volume. The frequency ujb and the amplitude Lb of the oscillation 
can be established from semiclassical arguments [a, @|. The former is usually 
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known as Bloch frequency and it is proportional to the applied electric field 
F, namely ub = eFa/h, where — e is the electron charge and a refers to the 
lattice period along the field direction. This periodic motion takes place in 
real and in k space, and it is known as Bloch oscillation (BO). Since scatter- 
ing processes, due to phonons or defects for instance, destroy the coherence 
necessary to support BOs, its experimental detection is highly nontrivial. 
The first BOs were detected in semiconducting superlattices 
More recently, cold atoms and Bose-Einstein condensates (BECs) in opti- 
cal lattices have been revealed as a very convenient scenario to observe BOs 
since scattering processes can be sig nificantly reduced Q H H- How- 
ever, decoherence effects cannot be removed completely and should be taken 
into account from a theoretical point of view. For instance, in BECs the 
atom-atom interaction gives rise to well-known dynamical instabilities which 
can destroy the coherence required to observe stable BOs I15|, lis ll7| • Still, 



there have been some experimental and theoretical proposals to avoid these 



undesirable effects [18|, [19| . Also, in organic molecules, which are very flexi- 
ble, the vibrations of the lattice are relevant and could rapidly degrade the 
electron quantum coherence. In this regard, the Peyrard-Bishop-Holstein 
(PBH) model of charge transport in DNA [20] was considered to demon- 
strate that polarons perform BOs even at realistic values of the carrier-lattice 



coupling [21 . 

In the last decade, a new interest in the dynamics of a quasiparticle 
affected not onl y by a DC field but also by a superimposed AC field has 



emerged [22|, |23|, |2J, |25| . Remarkably, it was demonstrated experimentally 
that a weakly interacting BEC in a harmonically driven tilted potential can 
support directed transport or large amplitude oscillations of the wave func- 

The latter are known as super 
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tion, depending on the driven frequency 
Bloch oscillations (SBOs) and their characterizing parameters can also be 
described semiclassically j27|. In this work we use the PBH model to study 
how the introduction of the carrier-lattice interaction affects the dynamics of 
polarons in this scenario. We will also calculate the current density associ- 
ated to the charge motion occurring under the applied fields and its Fourier 
spectrum to reveal the main frequencies involved in the polaron dynamics. 



2. Model 

The Hamiltonian of the PBH model can be written as 

H — %lat + %car + %int • 



(1) 
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The first term describes a one- dimensional anharmonic lattice according to 
the Peyrard-Bishop model |28j. A single degree of freedom X n IS cLSSl ffned to 
every site, taking into account its displacement from the equilibrium config- 
uration. %i at for a homogeneous lattice reads 



N 



(2) 



where m is the mass of each site and n = 1,2, ... N labels the sites along 
the system. Two potential terms appear in this Hamiltonian, namely a local 
Morse potential V(x n ) and a nonlinear anharmonic coupling between nearest- 
neighbors W(x n ,x n -i) 

,2 
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For the sake of concreteness, the fitting parameters are chosen according to 
those obtained to reproduce experimental DNA melting curves within the 
Peyrard-Bishop model, m = 300 amu, Vo = 0.04 eV, a = 4.45 A -1 , k = 
0.04 eV/A 2 and (3 = 0.35 A" 1 @. 

The charge carrier Hamiltonian 7{ car in the unbiased lattice is expressed 
within the nearest-neighbor approximation as follows 30 
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where T is the nearest-neighbor hopping and (V n and a n denotes the carrier 
creation and annihilation operators, respectively. Since there is not a com- 
mon value for the parameter T in the literature, we will take T = 0.1 eV 



hereafter as a good representative one |3l|, |32|, [33 



The last term in ([T]) takes into account a Holstein-like carrier-lattice in 



teraction as an on-site energy correction as follows [30 



n 
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(5) 



Here % denotes the carrier-lattice coupling constant. Ab-initio estimations 
of this coupling \ are scarce and therefore we will vary its magnitude in our 
numerical simulations. 
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In typical molecular systems it is possible to use a semiclassical approach 



due to the different time-scales of the charge and the lattice dynamics [34 
Thus, the dynamics of the carrier under a field with both DC and AC com- 
ponents such as F(t) = F + i^jsin(u;i + 4>) can be studied by way of the 
following Schrodinger equation 



. fc #n 

in—— 
dt 



= -[U + U x sin(wt + 0)]n^ n - T(ip n+1 + ^ n -i) + X^n , (6) 

where ip n is the probability amplitude for the charge carrier located at the 
nth site. The parameters C/j = eaFi with i — 0, 1 are the energy terms 
associated to the applied electric fields F i, a being the period of the lattice 
(e.g. a = 3.4 A in DNA). The last term in Eq. ([6]) describes the carrier- 
lattice coupling through the constant x an d the displacement x n from its 
equilibrium position. Newton's equations of motion for the displacements x n 
become 

d?x 

m —^ = -Vl^Xn) - W'(x n , x n -x) - W'(x n , x n+1 ) - x\^n\ 2 , (7) 

where the prime indicates differentiation with respect to x n . In what follows 
we will take the stationary polaron of the unbiased system under considera- 



tion 



30| as the initial condition for the integration of Eqs. ()6]) and ([7]) in a 



lattice subjected to the electric field F(t). 

3. Motion of the polaron in a biased system 



In Ref. [26| the dynamics of a BEC under superposed DC and AC forces 
was observed experimentally. Later, the authors of Ref. [27| studied the wave 
packet dynamics under these conditions in a rigid lattice (x = 0) by means 
of a semiclassical approach and numerical calculations. In both works it was 
reported that for a resonant AC field with frequency u = ug, the wave func- 
tion of the carrier performs a directed motion with oscillatory features with 
drift velocity depending on the phase <fi. On the contrary, if the frequency 
of the AC field is detuned with respect to w^, the wave packet exhibits a 
beating effect called SBOs j26[. These two situations will be referred to as 
resonant and detuned cases hereafter. In this work, we will focus on two set 
of parameters as representative values of these two different situations, for 
several values of the carrier-lattice coupling. 
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First, we analyze how the polaron evolves when the field F(t) = F + 
Fi sin(u;£ + 0) is applied. The time-domain evolution of the carrier wave 
packet obtained by direct integration of Eqs. (j6]) and (J7J) is shown in Fig. [1] 
for the resonant and the detuned cases, and x = 0.1 and 0.3eV/A. 



a) b) 




n 



Figure 1: Modulus of the carrier wave function in a lattice of N — 750 sites as a function 
of position and time with Fq = lOmeV/A F± — 0.5Fq and x = O.leV/A (left panels) 
and x = 0.3eV/A (right panels). The two upper panels correspond to resonant cases with 
uj = ujb and = 0. The two lower panels correspond to detuned cases with u = 1.2wb 
and 4> = 0- Light and dark regions indicate nonzero and zero values, respectively. 

Similarly to the case of the rigid lattice, the two upper panels of Fig. [1] 
(resonant cases with u = uib) show the directed motion of the polaron as 
well as superimposed BOs. On the contrary, the two lower panels of Fig. [T] 
(detuned cases with u = 1.2ub) shows that the carrier perform SBOs. It is 
to be noticed that by increasing the strength of the carrier-lattice coupling 
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the localization of the stationary states becomes larger and therefore the 
wave packet motion is more clearly observed (see Fig. [TU). We would like 
to stress that contrary to what happens in the rigid lattice (27| . due to the 
nonlinear interactions the wave packet remains localized in time except for 
long time-scales if the carrier-lattice coupling is large, as seen in Fig. [Tb). 

The dynamics of the polaron can be monitored in more detail by means 
of the centroid of the carrier wave function c(t) = n(t) — n(0) with n(t) = 
^2 n n\ip n (t)\ 2 . Similarly we define the dimensionless magnitude l(t) — £(t) — 
£(0) with = Yln=i n x nif) I a f° r the lattice displacements. We also study 
the Fourier transform of these magnitudes to reveal the frequencies involved 
in the dynamics of the polaron. 

Figure [2] shows that c(t) and l(t) perform a directed motion with short- 
time-scale oscillations in a lattice of iV = 2000 sites. The Fourier spectrum 
c(ui) indicates that the main frequency of such oscillations is equal to = 
eaF /h = 51.67 THz (very weak peaks at multiples of ojb are also found), 
while for l(u) the Morse frequency ojm = yZVoJm = 7.22 TH is the relevant 
one. Notice that the um is the harmonic frequency of the small amplitude 
oscillations of a mass m around the minimum of the Morse potential f l3a|) . 




t (ps) t (ps) 

Figure 2: Centroid of the carrier wave function (left) and l(t) (right) as a function of time 
in a lattice of N = 2000 sites for w = u>b, 4> — and \ — O.leV/A (upper panel) and 
X = 0.3eV/A (lower panel). Insets show the corresponding Fourier spectra. 

The results corresponding to the detuned cases (cj = 1.2ub) are shown in 
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Fig. EJ c(t) and l(t) display a more complex dynamics and absence of directed 
motion. In the case of c(t) the dynamics corresponds to the SBOs with a large 
period 2-k/Auj defined by the detuning frequency Au = 0.2ojb and a short 
period 2tt/ub- In addition, in c(u>) we obtained other peaks which can be 
described with the full analytical solution of the centroid motion in the rigid 
lattice [22I] [see the insets of the left panels of Fig. [3] . The main frequency 
involved in l(t) is again the Morse frequency but now we also observe the 
occurrence of a smaller peak at the detuning frequency [see the insets of the 
right panels of Fig. [3]. Therefore, we come to the conclusion that in both 
situations the oscillations of the lattice and the carrier are almost decoupled, 
at least at moderate applied fields. 




t (ps) t (ps) 

Figure 3: Same as in Fig.[2]but for detuned frequency (uj = \.2ujb and = 0) in a lattice 
of N — 750. 



4. Average current density 

In view of the oscillating behavior of the carrier wave packet described in 
the previous section, it seems reasonable to expect that the electric current 
behaves in a similar way. Therefore, and since the current is a macroscopic 
magnitude which can be directly observed in experiments, we calculate the 
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average current density J(t) according to the following expression 



20 



At) = ^J^J2 lm [r n ^n + i - Vn-i)] , (8) 

c n=l 

where m c is the mass of the carrier. 

As in the previous section, we will consider different values of the carrier- 
lattice coupling for the representative resonant and detuned case. Typical 
results of our simulations are collected in Figs. H] and [5j The left panels show 
the envelope of the average current density §E§ over a long time interval 
while the insets display the short time behavior for two different values of 
the coupling constant x — 0-1 eV/A and x — 0.3eV/A. 
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Figure 4: Left panels show the envelope of the average current density as a function of time 
in a lattice of N = 2000 sites for the resonant case for w = u)g, = and \ = 0.1 eV/A 
(upper panel) and \ — 0.3eV/A (lower panel). Insets show the short time behavior of 
the average current density. The right panels show the corresponding Fourier transform 
of the average current density. 

The average current density displays a well-defined oscillatory behavior 
in both cases, whose relevant frequencies match perfectly those obtained 
from the centroid motion in Figs. |2] and [3j Notice that these frequencies 
are in agreement to those analytically predicted for the rigid lattice in all 
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Figure 5: Same as in Fig.|4]but for detuned frequency (u — 1.2lub and = 0) in a lattice 
of TV = 750. 



cases (22J, |27| . We stress that the increase of the carrier-lattice coupling leads 
to a faster modulation of the average current density but, remarkably, the 
oscillations do not decay in time. 



5. Conclusions 

We have studied the dynamics of the carrier dynamics in the PBH model 
under superimposed DC and AC fields and found that it is similar to that 
found in the rigid lattice 



26, 27 



The carrier-lattice coupling leads to a 
distortion of the initial shape of the wave packet at long times and to a faster 
modulation of the centroid motion. Still, the polaron display SBOs which do 
not decay in time when the driven frequency is detuned with respect to the 
Bloch frequency. The carrier is partially decoupled from the lattice, whose 
relevant frequency is not the Bloch frequency but the one associated to the 
Morse potential. 
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